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Abstract 

The results of a recent experiment measuring the reaction pp dK~^K^ near threshold are 
interpreted in terms of a spectator model that encapsulates the main features of the observed 
K^K^ invariant mass distribution. A fit to this data leads to an imaginary part of the isovector 
scattering length in the KK channel of Im{ai) = —(0.63 it 0.24) fm. We then use the Flatte 
representation of the scattering amplitude to infer a value Re{ai) = —(0.02 it 0.02) fm for the 
real part under the assumption that scaling (Baru et a^li^l) is approximately satisfied. We show 
further that it is not possible to exclude the eff'ects of ■K'^rj to K'^K^ channel coupling within the 
context of our model. 
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1. Introduction 



The reaction pp — > dK^K^ has been measured at an excess energy of Q = 46 MeV with 
the spectrometer ANKE at the cooler synchrotron COSY-Jiihch[l]. This experiment has 
achieved a resolution of 1 — 8 MeV for the invariant mass-spectrum of the KK system which 
considerably improves the data base. Therefore it is timely to investigate what information 
regarding the KK scattering length can be extracted from this data. 

A model-independent analysis of the angular distributions shows a dominance of an s wave 
between the two kaons accompanied by a p wave deuteron with respect to the mesons. This 
has been interpreted as evidence for a dominant production of the ao(980) near threshold 
P, 1^ . Detailed calculations of the reaction pp dK^K^ performed by Grishina et al. j3| 
confirm this conclusion. In particular the non-resonant K~^K^ production is shown to be 
small below excess energies of Q = 100 MeV. Fadeev calculations for the K~d reaction have 
extracted a scattering length ax-d = (1-80 — il.55) fm (using the Bethe sign convention 
here)|4l. A summary of K~d scattering lengths obtained by different authors can be found 
in Ref. [5| with Re{aK-d) ranging from 0.01 to 1.92 fm. 

Data obtained more recently at COSY have been analysed by Sibirtsev et al. Q. These 
authors have extracted limits for the Kd scattering length from the data of Ref. ^ and 
found the magnitudes of both the real and imaginary parts of this quantity to be less than 
1.3 fm. A chiral approach to the reaction pp — > dK^K^ has been formulated by Oset, Oiler 
and Meissner j|6|] that suggests a potentially large sensitivity of the data to the Kd scattering 
length. A study [?[ based on this formalism finds an optimal fit to the COSY data for a 
purely imaginary Kd scattering length of = (0.0 + il.O) fm. 

The finding of a small real part of the scattering length aj^^ in Refs. js], 0| encourages us 
to explore the working assumption that the deuteron's role is simply that of a spectator, 
and that the K^K^ production in pp — > dK^K^ is dominated by final state interactions in 
the meson-meson channels. Since the reaction pp d-K^i] has a lower threshold than pp — >■ 
dK~^K^, both 7r+?7 production and K^K^ production have to be treated simultaneously. 
We therefore consider a coupled channel description of the two-meson final states t^^t] 
and K^K^. In the vicinity of the K^K^ threshold, one expects the T^^rj production to 
dominate. This allows one to make a simple estimate of the two kaon production via the chain 
Ti^T] ^ ^ K^K^ . In order to include the other reaction chain K^K^ — > ag" — > K^K^, 
one would need a model for the production operators. In view of this we assume that the 
first reaction chain is dominant, but explore the influence of the second chain empirically. 
This approach cannot make any predictions about the three additional measured angular 
distributions or the dK^ invariant mass distribution reported in Ref.jH. 



2. Spectator model calculation 

The scattering amplitudes for K^K^ — ^ K^K^ and tt^t] — > K^K^ are related by two- 
channel unitarity. Therefore at low relative momentum k in the two kaon channel both 
amplitudes can be expressed in terms of the KK isovector complex scattering length ai = 
as — iPs- Calling these two channels 1 and 2 respectively, one has |^ 

Tn = r] e^'^^ - 1 ^ -2ikai{l - ikai) 

e-^-Ti2 = t^l-v' ^ 2i^{l - I3sk) (1) 
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where 61^2 are the phase shifts in the corresponding channels, Su = ^1 + ^2, and rj is the 
inelasticity. 

We now assume that the s wave single differential cross section da^'^^ / cLPq for the pp — >■ 
dK^K three-body decay can be written as the proportionality 

da^'^dPo ~ r™'' dt\AT^2 + BTu\^ ^ Csk\{l - P,k) + pe'^^P^k pAt (2) 

where p exp{i(f)) = iB/A{l + ias/ j3s) exp{—iSi2) summarizes the admixture of channel 1 
via the ratio of the production operators A and B for channels 1 and 2, At is equal to the 
difference in limits of the phase space integration, 



At = t„iax - tmin ~ \/ {s - {niD + i^o)^)(s - {mo - -Pq)^), (3) 

-s = {pa + PbY is the invariant mass of the colliding protons, and Cg a coefficient of propor- 
tionality. In principle both Cg and p expiicl)) will depend on specific details and kinematics 
of the production mechanism in the two channels in question. We assume these factors to 
be slowly-varying for simplicity. For the three-body collision the KK invariant mass is 
restricted to the range from threshold (Po)mm = {Mk+ + Mf^o) to {Po)max = (v^ — triD) so 
that At vanishes at the upper limit; m^i is the deuteron mass. At the COSY proton beam 
energy of Tp = 2.65 GeV {pp = 3.46 GeV/c) at which the experiment was performed, these 
limits are (Po)mm = 0.9914 GeV and iPo)max = 1-038 GeV resp. This restricts the K+K^ 
and ir'^T] CM momenta. 



k = -Po 
2 " 

1 



{Mk+ - Mjiofl 1/2 r {Mk++Mko) 



p2 



p2 



1/2 



P 







P=-^Pol- ' l- ^--^-' (4) 







1/2 



to lie in the very limited intervals < A; < 0.15 GeV and 0.33 < p < 0.36 GeV respectively. 
The Po dependence of the cross section is thus completely determined by the variation of k 
and At with Pq as fixed by Eqs. (jD) and Q. The scale constant Cs[ph GeV~^] controls the 
absolute value of the s-wave cross section while the parameter f^g = —Im{ai) determines 
its skewness in the case of no admixture, p = 0. 

The expression equivalent to Eq. Q for p waves to leading order in k is 

da^P^/dPo ^ Cpk^At (5) 

which for Cp = 150.6pb GeV~^ reproduces the corresponding contribution determined in 

n 

The fits to the various data shown in Ref.jl| employ a six-parameter transition matrix 
fitted simultaneously to the two invariant mass distributions K'^K'^ and dK^, plus 
three angular distributions. In order to compare the cross section given by the spectator 
model with their result for the K~^K^ invariant mass distribution, we have determined the 
pair of s wave constants [Cs,/?s] by performing a fit to the data points in Fig. [T] that lie 
above threshold, using the sum of Eqs. (j2I) and (jS)) after expanding the former up to order 
0{k), initially for the case of no mixing, p = 0. In Sect. 4 we discuss the effect of including 
channel admixtures as well as relaxing the assumptions on Cg. 
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One finds Cs = 22.03/i6 GeV~^, together with an imaginary part of the scattering length 
of = Im{—ai) = 0.63 fm at Xmin — 0.65. The fit to experiment for this parameter set 
is shown in Fig. 1. By integrating out over the allowed range of Pq, one obtains a total 
cross section of a{pp dK^K^) = 36 nb, in agreement with the experimental cross section 
[38 ± 2{stat) ± U{syst)] nb reported in Q. 



3. Flatte parametrization for scattering length 

The spectator model assumption Eq. (0) can only fix the imaginary part of the scattering 
length for given [p, 0] : it does not allow one to extract the real part directly. In order to 
make further progress one therefore has to introduce specific model dependent assumptions 
regarding the structure of the KK scattering amplitude in the isovector channel that describe 
the formation and decay of the ao(980). 

There is little experimental information available on the KK scattering length. However, 
within the last decade the production of ttt] pairs have been studied experimentally in proton- 
antiproton annihilation, pion-proton reactions and in $ decay, and analysed using either 
Flatte distributions jnl, [ij, 3, or other approaches [isl. Ea Ht^ . The Flatte scattering 



amplitude ^ISii] near the KK threshold takes on the form 

c - ~ 1 ^dKRk , . 

-eao + 2i9KKk + 9nvP) 

where k and p are the CM momenta defined in Eqs. (0)), while the ao(980) partial decay 
widths = (Jt^^p and Tj^j^ = QkR^ define the dimensionless coupling constants {g-nr^^ Qkk)') 
= Mao ~ + Mj^o) is the (real or virtual) binding energy of the ao(980) relative 

to the KK threshold. Eq. © leads to the following result for the KK isovector scattering 
length 111 

\9KK R -1 fr7\ 



when rewritten in terms of the dimensionless ratios |20| R = gKR/Qirri and a = 2eag/TT,ri 
instead of the individual coupling constants and energy parameters. Here pth = p[{Po)min] = 
0.33 GeV = 1.65 fm~^ is the threshold momentum for KK production in the ttt] channel. 

The R and a are seen to be invariant under a scahng gj^j^ —> Xg^R ^tc. of the original 
Flatte parameters. Table I shows a compendium of mass, width and coupling constant 
values for an (980) taken from the recent literature 11, l^, IE 3 3 3, 3 ^^"^ partiall 



summarised in jl9|, |20| . We have also used Eq. (|7j) to extend the calculations given in 
for the scattering lengths associated with these paramaters. The results are displayed in 
order of increasing a in Table I. 

From Eq. ((Tj) one sees that there is a linear relation between the real and imaginary parts 
of the (Flatte) scattering lengths, 

Re{—ai) = Im{—ai)a (8) 

If we now take Im{—ai) = 0.63 fm from the spectator model fit, this fixes the slope of 
the straight line in Eq. (jH)). The result can then be compared with the values of Re{—ai) 
calculated from Eq. ^ at the values of a and R given in Table I. This comparison is 
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shown in Fig. 2. The spectator model result is seen to be consistent with all of the data 
(except perhaps that of Ref. once the uncertainties in a and hence Re{—ai) due to the 
experimental error in the an (980) mass are taken into consideration. 

It has also been pointed out[20| that direct measurements of the KK elastic scattering 
cross section near threshold would serve to determine both R and a independently, and 
thus the scattering length as given by Eq. ((7j). Such data is not available yet. However, by 
assuming that scale invariance is valid, one can infer a value for a to be used in Eq. (jS)) 
by arguing as follows: the combinations a in Table I would have appeared as a common 
constant etc had the data satisfied scaling exactly. It is therefore reasonable to estimate this 
unknown constant by the mean of the a's given in Table I, weighted by the inverse of the 
average error introduced into each a by the uncertainty in the ao(980) mass value used to 
calculate it. The resultant mean and its standard deviation is then = 0.031 ± 0.054. By 
also considering the changes in slope of the linear relation Eq.(jH)) due to the errors introduced 
into Re{—ai) by the mass value uncertainty one obtains the order of magnitude estimate 
f5s = (0.63 ± 0.24) fm for the error in (3s. Inserting these values for j3s and Oc together with 
their estimated errors into Eq. (jH)) gives 

ai ^ -[(0.02 ± 0.02) + i(0.63 ± 0.24)]fm (9) 

for the isovector scattering length. 

Note that this result relies on two very different sources of information: (i) a direct 
determination of /m(ai) based on the spectator model fit of the pp — >■ dK^K^ data for p = 0, 
and (ii), an inferred estimate of Re{ai) based on the assumption that scaling is approximately 
valid for the experimental parameters extracted assuming a Platte distribution. 

4. Effect of mixing with K+K^ K+K^ 

We now repeat the fit of the data points in Fig. 1 in the presence of mixing. Since 
p > 1 can lead to unphysical values of f3s < 0, we illustrate the effects of mixing for the two 
limiting cases labelled [p, 0] = [1, tt] and [1, 0] in Fig. 2. The corresponding values of (3s swing 
from 0.21 to 7.6 fm respectively as illustrated by the gradients of these two straight lines. 
However, by choosing an intermediate value of appropriately in the interval < < tt for 
given p, in this case [p, 0] = [1, 110° ± 8°], one again recovers a fit (Xmm — 0.82) to the data 
in Fig.l at the indicated value of (3s given in Eq. Q that is essentially indistinguishable 
from the no-mixing case. Similar observations hold for any < p < 1 except that the swing 
in possible values of (3s is smaller. 

Another source of uncertainty is associated with the lack of a detailed description of the 
production mechanism itself. As an illustration of this we compare with the specific model for 
the production mechanism suggested in Ref. [10] that codifies selection rules explicitly. One 
finds that their scattering amplitude, which also considers resonance production, reduces to 
the first term of Eq. (j2)) at low momentum with Cs — > (CsP'jj), where pi? is the deuteron 
momentum. Making this replacement in Eq. (0), one obtains an equally good fit (Xmm ~ 
0.20) of the mass distribution in Fig.l with (3s = 0.63 fm by changing the mixing angle from 
110° to 59° ±4°. 

At the level of the spectator model as envisaged in this note it is therefore impossible 
to decide whether or not channel mixing is important without introducing model specific 
assumptions for the production operators A and B. 
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Invariant mass (GeV) 



FIG. 1: The invariant mass cross section for the production of K^K in pp dKK. The soHd 
hne denotes a fit to the dataj3] using the spectator model without mixing (the fits in the 
presence of mixing are essentially identical for reasonable choices of the mixing parameters [p, (p] ) . 
The broken curve shows the corresponding fit calculated in 0] using a six-parameter transition 
amplitude squared. Shown separately are the p wave contributions to the cross section obtained 
from Eq. © (solid curve) , and the fit (broken curve) given in P| . The error bars show the statistical 
uncertainties of the data. 
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FIG. 2: Plot of Re{—ai) (filled circles) versus the dimensionless parameter a = 2eo/TT^f^. The 
plotted points have been tagged by the corresponding references in Table I. The horizontal and 
vertical error bars show the errors introduced into a, and hence Re{—ai), by the uncertainties in 
the associated ao(980) mass values. The heavy straight line shows the calculated value of Re{—ai) 
using Eq. © with j3s = Im{—ai) = 0.63 fm. The label [1, llO'^] gives the mixing parameters that 
lead to the same value of (3s upon including channel mixing as characterised by the parameters 
[p, 0] defined in the text, while the remaining two broken straight lines labelled [l,vr] and [1,0] 
show the minimum and maximum values that (3s can assume for p = 1- 



9 



TABLE I: Summary of calculated values of the isovector scattering length oi taken from Eq. ((Tj), 
arranged in order of increas ing a = 2eo/TTjj^, for the various mass, width and coupling constants 
for the ao(980) as extracted ,2u] from the literature. 



Rp f 


M (MpV) 


F (MpY) 


r ('MeV'l 




yKK 


R 




n -1 






984.8l:|:^ 


-6.6 


121.5 


0.373 


0.412 


1.10 


-0.108 


+0.071 - 


0.66? 


[12] 


992+1^ ' 


0.6 


145.3 


0.446 


0.560 


1.252 


0.008 


-0.006 - 


0.76i 


M 


995tll 


3.6 


125 


0.384 


1.414 


3.685 


0.058 


-0.13 - 


2.23i 




9991^ 


7.6 


143 


0.439 


0.516 


1.176 


0.106 


-0.075 - 


0.70i 


[12j 


W03tll 


11.6 


153 


0.470 


0.834 


1.776 


0.152 


-0.16 - 


1.05i 


[17] 


999^1 


7.6 


69 


0.212 


0.222 


1.048 


0.220 


-0.13 - 


0.61i 


[13] 


lOOll^ 


9.6 


70 


0.215 


0.224 


1.043 


0.274 


-0.16 - 


0.59i 


m 


998l{0 


6.6 


46 


0.141 


0.476 


3.376 


0.287 


-0.54 - 


1.89i 



10 



